Conductivity of a quasiperiodic system in two and three dimensions 
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A generalization of the Aubry- Andre model in two and three dimensions is introduced which allows 
for quasiperiodic hopping terms in addition to the quasiperiodic site potentials. This corresponds 
to an array of interstitial impurities within the periodic host crystal. The resulting model is exactly 
solvable and I compute the density of states and the ac-conductivity. There is no mobility edge as 
in completely disordered systems but the regular ac-conductivity and the strongly reduced Drude 
weight indicate a precursor of the Anderson transition as the Fermi energy goes from the center to 
the band edges. 
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The description of electron motion in a non-periodic 
potential within the single-particle approximation is an 
elementary problem in solid state physics and still un- 
solved for the case of a macroscopic number of impurities 
which break discrete translational invariance. Especially 
the phenomenon of Anderson localization in a tight- 
binding model with uncorrelated random site potentials 
defies exact analytical treatment up to now. Therefore, 
several attempts have been undertaken to study mod- 
els with quasiperiodic potential, which constitute a case 
intermediate between the perfect crystal and a fully dis- 
ordered system. The simplest realization of this situation 
is a tight-binding hamiltonian Hi where the site energies 
vary in a quasiperiodic fashion: 



Hi = H Q + ^ 



|n)(n| 



(1) 



The |n) are local orbitals at site n of a d— dimensional 
hypercubic lattice, Hq describes the hamiltonian of the 
periodic lattice and g m is a quasiperiodic function of m. 
A possible choice for g m is 



fjr, 



2g cos(Q • m) 



(2) 



with a vector Q whose components Qj are incommensu- 
rable with 7r. This is the d— dimensional Aubry-Andre 
model and unsolvable for d > 1 despite its simplicity. 
Even in d = 1 only very limited results can be obtained 
exactly, e.g. the localization property of the eigenstates 
if 5 > 1-- 

A general criticism applies to all models of type |T]) : The 
doping of the periodic host lattice with impurities is usu- 
ally not of the substitutional type but the additional 
atoms sit at interstitial positions. They modify not only 
the site energies but change locally the orbitals and in 
turn the overlap integrals defining the hopping matrix 
elements. Furthermore, an empty impurity orbital may 
serve as intermediate state for a hopping process connect- 
ing distant sites of the host lattice. These effects may be 



phenomenologically taken into account by the following 
generalization of fl} : 



H = Hi+J2 ^^(n-m,n'-m)|n)(n'| 



(3) 



The function g' m (n — m, n' — m) depends on the position 
of the impurity m and decays with growing distance of 
the sites n, n' from m. This generalization of {1} seems 
to be more complicated than the Aubry-Andre model if 
both the g m and the g' m vary quasiperiodically with m. 
It is therefore surprising that one may construct such a 
model which is completely analytically solvable and has 
a solution in arbitrary dimensions. "Complete solution" 
means here that the hamiltonian is diagonalizable and 
computation of the transport properties can be done an- 
alytically. As in most non-trivial soluble models, the solv- 
ability rests here on a relation between the quasiperiodic 
functions g m and g' m . The starting point is the periodic 
hamiltonian 



H 



E 

(n,n') 



|n)(n'|+h.c. 



(4) 



where (n, n') denotes a pair of next neighbors on the 
d— dimensional hypercubic lattice A = Z d with lattice 
constant a = 1. Its matrix elements read in momentum 
space 



H (p,p' 



-2t^2cos( Pj )S(p 



(5) 



with p, p' G [— 7r, 7r] d . 

We introduce two types of potentials related to the site 
m: 

I. A local potential term with matrix element V m (n,n') 
in position space: 



V m (n, n') = g m (2ir) d ]J S( nj - n' 3 )S( nj - to,). (6) 
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II. A non-local term: 



O \ d d 

l) n 



C(n,n') = 



(-1)' 



L ( % .- mi -l/ 2 )( n ;.-m, -l/2) 



(7) 



The functional form of ([7]) agrees with the general ex- 
pression in and can be interpreted as perturbative 
description of a hopping process, where the electron goes 
from site n first to the impurity site m + 1/2 on the dual 
lattice and from there to the site n'. The technical rea- 
son for this choice of the non-local term is the fact that 
now both potentials have the same form in momentum 
space: 

Kn(p,p') = ff m exp(im-(p'-p)) (8) 
C(P,P') = <4exp(ifii-(p'-p)). (9) 

with rhj — rrij + 1/2. Because both forms of the potential 
factorize with respect to the dimension index j, also a 
mixed type is possible, which is for some j of type I and 
for the rest of type II. All of these possibilities can be 
parameterized through a vector m with integer and/or 
half-integer components, i.e. m 6 A' = (Z/2) d and we 
may drop the primes in ([9]). The generalization of the 
Aubry- Andre model in d dimensions reads then 



H = H Q + 



(2Tr) d ^ 



and 



2gioc cos(Q • m), meA 
2g n ioc cos(Q • m), m G A' \ A 



(10) 



(11) 



Here, Q denotes some vector in M. d such that Qj/ir are 
irrational. We may further assume that the numbers 
Qj/ft are linearly independent over Z. In this case, no 
point of the lattice is equivalent to another. The original 
Aubry- Andre model is recovered for g n \ oc = 0. We set 
in the following g\ oc = <? n i oc = .<?■ The potential term 



space: 



V m can be rewritten in momentum 



Vq(p, p') = g (Mp' - p + Q) + Mp' - p - Q)) • 

Here, 5^{x) denotes the delta-function with fundamen- 
tal period 4-7r: 8&k(x + Ann) = 5^(x) for n G Z. The 
quasiperiodic potential translates in momentum space 
to displaced delta-functions, which however do not have 
2tt— periodicity but 47r— periodicity due to the addition of 
the non-local terms of type II. As a result, only finitely 
many points in momentum space are connected through 
the potential term, even if Q/tt is irrational. 
The simplest case is realized if just two sites in momen- 
tum space are connected. This happens if one of the Qj 
lies in the interval [tt, 2tt\. The other components of Q 
can be chosen at will. The Brillouin zone splits into three 



regions Rq,Ri and i?2, where Rq contains points, which 
are not affected by the potential (fT2")l at all, whereas Ri 
and i?2 are mutually connected through the potential: 



Ri = 



Ri = 



Rn = 



i[[-*,*-Qi] 

7=1 
d 

IjHr + Qj.Tr] 

7=1 
d 

l[{-ir,ir}\(R 1 UR 2 ) 



(13) 



(see Fig Q] for the case d = 2). 

It has to be emphasized that this splitting is equivalent 




FIG. 1: The three regions Rq, R lt R 2 for d = 2. 

to a reduction of the Brillouin zone and concomitant gen- 
eration of sub-bands only if Q/7r is a rational vector. In 
the general case the model does not have a band struc- 
ture. 

Nevertheless, vectors p from Ri can be used to label 
the two-dimensional invariant subspaces TL P , which are 
spanned by |p) and |p + Q). On these subspaces, the 
hamiltonian can be easily diagonalized: 



H 



-2tJ2j C0S (Pj) 9 

H p \ 9 -2<E J C0S (Pj - + Qj) 

(14) 

In the following we have set for simplicity t = 1/2. Fig. [5] 
and [3] give the density of states in two and three dimen- 
sions, respectively. The total DoS is composed from two 
"bands", where the p— band stems from the region Rq, 
which is unaffected by the quasiperiodic potential and 
the eigenstates of (1101) are the pure momentum eigen- 
states |p) for p e i?o- These are located in a region 
around E = 0, the band center. The q— band is gapped 
and contains the eigenstates \tpp) in H p for each p G i?i 
and energy eigenvalues Ep. It is located predominantly 
at the band edges. 
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becomes periodic in fc— direction. FigfJ] shows D%/f3 for 
small temperatures and Ep located in the gap of the 
q— band as function of Q x and Q y . The Drude weight 



kTD 




FIG. 2: Density of states for d — 2. The potential is quasiperi- 
odic with Q x — 2tv — 4 and Q y = 2-k — 2\/2. The coupling 
constant is g — 2. 
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FIG. 3: DoS for d = 3. Here, Q x = 2-k - 4, Q y = 2tt - 2^3 
and Q z = 2-k — \/E with coupling g — 2. 



To compute the conductivity, we first note that the 
p— band is not affected by the potential: The current 
commutes with the hamiltonian and we get a diagonal 
ac-conductivity (f3 = (fcT) -1 ): 

Re(a kk (oj))=TrD (f3)S(u J )+a r k l 9 (oj) (15) 

with vanishing regular part, a^ 9 (u>) = and the Drude 
weight Dq(/3) of the pure system, which is only reduced 
because Rq does not cover the whole Brillouin zone. 
In the q— band, however, we obtain a non-vanishing <tI? 9 
because the current J does not commute with H on the 
spaces Ti v : 



Jk 



sin(p fe ) 

sin(p fc + Q k ) 



(16) 



The (anisotropic) Drude weight D%(j3) of the q— band is 
computed as 



(17) 



Dl(f3)=p dp n F (E±)\(^\J k \^)\ 



Ri 



with the Fermi- function np{E) = (1 + exp /3(E — Ep)) 1 . 
Dt approaches Dq for Q k — > 0, 2ir, i.e. when the potential 



FIG. 4: Drude weight of the g-band D%/@ for d = 2. The 
maximum at Q a = 7r and = corresponds to the Drude 
weight of the pure system Do- The minimum along Q y = 2tv 
is due to the vanishing area of Ri in this case. 

D? approaches a nonzero large coupling limit q — > oo 
with the exception of the (rational) value Qk/^ = 1, 
where the expectation value of Jk in the eigenstates of H 
goes to zero. This behavior is seen in Figj5] Apart from 





[31,31,31] 




[2.28,2.81,4.05] 





FIG. 5: Drude weight D%//3 for d = 3 and two different 
vectors Q. The large coupling limit is zero for periodic 
Q = (n, n, 7r) and nonzero for quasiperiodic Q = (2tt — 4, 2tt — 
2V2, 2-it-VE). 

D^, which is reduced in the q— band due to scattering 
from the quasiperiodic potential, the regular part of the 
ac-conductivity does not vanish: 



req / \ 



1 - e- f3uJ 



dp n F {E-)\^~\ J k \iPZ)\ 2 S(Lo + E~- £+; 



(18) 



E p is bounded from below 



with Ep < 0. Because E£ 

by 2g, the ac-conductivity vanishes for frequencies below 
this threshold. a^% a (uj) is plotted in Fig [6] for various 
values of /?, g and Q in d = 2. The form of a re9 (ui) in 
three dimensions is very similar. 

We have introduced and solved a tight-binding model 
with quasiperiodic local and non-local potentials in ar- 
bitrary dimension. The eigenstates fall into two groups, 
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FIG. 6: a^% 3 (u) for various parameters: C3 corresponds to 
g — 2,(3 = 1,Q = (2tv — 4, 2tv — 3), C4 has the same g and 
Q but (3 = 5, C2 and Ci both have g = 2.5 and (3 = 1 but 
different Q: (2tt-4, 2tt-3) and (2tt-4.8, 2tt-2.5). Although 
Tli has the same area in both cases, the ac-conductivity is 
quite different. Ef = — 2 in all cases; the Fermi energy lies 
within the lower q— band. 

where one group (states in the p— band around E = 0) 
exhibits the properties of the pure system but states 
at the band edges ((/—band), although still extended, 
have different transport properties due to the influence 
of the quasiperiodic potential: a strongly reduced Drude 
weight and non-zero regular ac-conductivity. This can 
be interpreted as a precursor of the Anderson transi- 
tion expected for lattices with uncorrelated site poten- 
tials. For irrational values of Qj/tt the system has no 
band structure in the usual sense. Nevertheless, the den- 
sity of states is absolute continuous and therefore "band- 
like"—. This feature has been found numerically in the 
two-dimensional labyrinth-tiling^, the analogous three- 
dimensional model^ and the square Fibonacci tiling^ for 
a certain range of parameters. The "band-like" spec- 
tra are prima facie closer to the experimentally observed 
quasiperiodic systems' than the models which exhibit 
a singular continuous spectrum like most of the one- 



dimensional examples 6 . But this does not mean that the 
present model is typical for real quasicrystals. The main 
objection to the physical relevance of the form for 
the non-local term is the following: the hopping matrix 
elements decay algebraically with distance from the im- 
purity site whereas the decay should be exponential, fn 
this respect the model resembles the Lloyd model, which 
describes uncorrelated disorder with a broad (lorentzian) 
distribution of the site potentials 7 . However, only the 
disorder-averaged DoS can be analytically calculated in 
the Lloyd model, not the transport properties. It is there- 
fore only partially solvable. Of course, the "long-range" 
nature of the quasiperiodic hopping term has a delocaliz- 
ing effect on the eigenstates and is probably the reason for 
the absence of a mobility edge in the generalized Aubry- 
Andre model. 

Besides the conductivity, one may study the temporal 
behavior of wave-packets initially located at a single site 
and the associated (anomalous) diffusion of the electrons. 
This has been done numerically for several quasiperi- 
odic systems in two and three dimensions'". Moreover, 
some theoretical predictions relating spectral and diffu- 
sive properties have been made'. These predictions can 
now be tested analytically in the present model, which 
will be the subject of a forthcoming paper. A second line 
of future investigation is a perturbation theory around 
the exactly solvable point g\ oc = g n i oc - This introduces a 
2ir— periodic displaced delta- function with a weight pro- 
portional to <7i oc — gnioc, a quantity which may serve as 
small parameter of the perturbative expansion. It is pos- 
sible that localized states appear as soon as gi oc — ffnioc be- 
comes nonzero, in which case the localized regime would 
be perturbatively accessible. 
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